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Abstract 

We discuss the general framework for the construction of new models containing 
a single, fermion number zero scalar leptoquark of mass ~ 200 — 220 GeV which can 
both satisfy the DO/CDF search constraints as well as low energy data, and can lead 
to both neutral and charged current-like final states at HERA. The class of models 
of this kind necessarily contain new vector-like fermions with masses at the TeV scale 
which mix with those of the Standard Model after symmetry breaking. In this paper 
we classify all models of this type and examine their phenomenological implications 
as well as their potential embedding into SUSY and non-SUSY GUT scenarios. The 
general coupling parameter space allowed by low energy as well as collider data for 
these models is described and requires no fine-tuning of the parameters. 
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1 Introduction and Overview 



1.1 Current Status of the Leptoquark Scenario 

The observation of a possible excess of neutral current (NC) events in e~^p collisions at 
high— by both the H1[|T| and ZEUS0 Collaborations have sparked much fervor in both 
the theoretical and experimental communities. This excitement has now been heightened 
by the recent announcement that both experiments may also be observing a corresponding 
excess in the charged-current(CC) channelQ. If these events are not merely a statistical 
fluctuation, it is clear that new physics must be invoked in order to provide a suitable expla- 
nation, e.g., compositeness appearing in the form of higher dimensional operators[Q, exotic 
modifications of the parton densities^, or the resonant production of a new particle0, 
such as a leptoquark (LQ) or squark in supersymmetric models with R-parity violation. 

If the excess is resonant in the x distribution 0, a popular interpretation ||^, |^ invoked 
in the NC case is the s— channel production of a ~ 200 — 220 GeV scalar {i.e., spin-0) 
leptoquark with fermion number {F) equal to zero. These quantum numbers arise from the 
requirements that (z) the observed excess appears in the e^p rather than the e~p channel, 
[a) the Tevatron search constraints^ exclude vector (spin-1) leptoquarks with masses near 
200 GeV, and {in) any discussion of leptoquark models has been historically based on the 
classic work by Buchmiiller, Riickl and Wyler (BRW) 0] . In that paper the authors provide 
a set of assumptions under which consistent leptoquark models can be constructed; these we 
now state in a somewhat stronger form: 

(a) LQ couplings must be invariant with respect to the Standard Model (SM) gauge 
interactions, 

(b) LQ interactions must be renormalizable, 

(c) LQs couple to only a single generation of SM fermions, 

(d) LQ couplings to fermions are chiral, 
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(e) LQ couplings separately conserve Baryon and Lepton numbers, 

(f) LQs only couple to the SM fermions and gauge bosons. 



Amongst these assumptions, both (a) and (b) are considered sacrosanct whereas (c)-(e) are 
data driven by a host of low energy processes. Assumption (f ) effectively requires that the 
leptoquark be the only new component added to the SM particle spectrum which seems quite 



unlikely in any realistic model. Based on these classical assumptions it is easy to show||10|| 
that all F = scalar leptoquarks must have a unit branching fraction into a charged lepton 
plus jet {i.e., Bi = 1). This lack of flexibility presents a new problem for the leptoquark 
interpretation of the HERA events for two reasons: (i) leptoquarks with = 1 clearly 
cannot accommodate any excess of events in the CC channel at HERA since these would 



require a sizeable leptoquark decay rate into neutrino plus jet, (ii) both CDF|]I2| and D0p3[| 
have recently presented new limits for the production of scalar leptoquarks at the Tevatron 
using the next-to-leading order cross section formulae of Kramer et al. [T^ . In particular, in 
the eejj channel, DO finds a 95% CL lower limit on the mass of a -B^ = 1 first generation 
scalar leptoquark of 225 GeV. DO has also performed a combined search for first generation 
leptoquarks by using the eejj, evjj and vvjj channels. For fixed values of the leptoquark 
mass below 225 GeV, these search constraints can be used to place an upper limit on 5^. 
For MiQ=200(210,220) GeV, DO obtains the constraints B^ < 0.40(0.62,0.84) at 95% CL. 
Of course if CDF and DO combine their searches in the future, then the 225 GeV bound may 
rise to ~ 240 GeV, in which case even stronger upper bounds on Bi will be obtained. 

Besides the obvious need to provide an potential explanation for the HERA data which 
satisfies all other experimental constraints, it is perhaps even more important to explore in 
a more general fashion how one can go beyond the rather restrictive BRW scenarios. Even 
if the HERA events turn out to be statistical fluctuations, we will show that by the removal 
of the least tenable of the BRW assumptions, we can find important ways to extend the 



possible set of leptoquarks that may be realised in nature. Since, as was mentioned above, 
it is difficult to believe that the addition of the leptoquark would be the only extension to 
the SM spectrum in any realistic model containing such a field, it is clear that assumption 
(f) should be abandoned. We now explore the consequences of this possibihty. 

1.2 Enlarging the Framework of Leptoquark Models 

In order to satisfy all the experimental constraints it is clear that we need to have an F = 
scalar leptoquark as before, but now with a coupling to SM fermions given by 

jO-wanted = [^^u^u" + Xded""] ■ LQ + h.c. , (1) 

with comparable values of the Yukawa couplings and A^. This fixes the Icptoquark's 
electric charge to he Qlq = ±2/3; no other charge assignment will allow the leptoquark to 
simultaneously couple to ej and uj as is required by the combination of HERA and Tevatron 
data. An alternative possibility, if neutrinos are Dirac particles, or if i/^ is fight and appears 
as missing px in a HERA or Tevatron detector, is the interaction 

^Lnted = [K^'u + X'.e'^d] ■ LQ' + h.c. (2) 

It is important for later analysis to note that these two interactions cannot simultaneously 
exist as the BRW assumption (d) above would then be strongly violated. Unfortunately, both 
of these Lagrangians as they stand violate assumption (a) above, in that they are not gauge 
invariant with respect to SU (2)^. We must then arrive at one of these effective interactions 
indirectly by some other means than by direct fundamental couplings. In order to do so it 
is clear that we must be willing to abandon at least one of the BRW assumptions (a)-(f) 
and it is evident that (f) is the one most easily dismissed. Hence we will assume that the 
leptoquark has additional interactions besides those associated with SM gauge interactions 
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and the Yukawa couplings to the SM fermions. We note, however, that fine-tuning solutions 
can be found which allow the assumption (c) to be dropped as a condition that applies in 
the mass eigenstate basis; these will not be discussed in detail here although it is important 
to understand how fiavor mixing plays a role in leptoquark dynamics in realistic models. 

In principle there are several alternatives as to what kinds of new additional inter- 
actions one can introduce, two of which we now briefiy discuss. In a recent paper, Babu, 



Kolda and March- Russell [jT5| considered an interesting model with two different leptoquark 
doublets, one coupling to Ldf^ and the other to Lu'^ (with L being the SM lepton doublet). In 
this model the electric charge Q = 2/3 members are mixed through a renormalizable coupling 
to the SM Higgs field with the mixed leptoquarks forming mass eigenstates that can couple 
to both ej and vj as desired with the ratio of strengths controlled by the amount of mixing. 
The new interactions in this case are quite complex and a certain amount of fine tuning is 
necessary to get the spectrum and couplings to come out as desired. The rich phenomenol- 
ogy of this scenario, which now involves four leptoquark mass eigenstates of various charges, 
should be further studied in detail. A second scenario has only been briefiy mentioned in 
the recent paper by Altarelli, Giudice and Mangano[|l6| who considered the possibility of 



at least temporarily violating both conditions (a) and (b) via non-renormalizable operators. 
These authors show, however, that both (a) and (b) can be restored by the introduction of 
new heavy fermions to which the leptoquarks couple in a gauge invariant fashion and which 
are then integrated out to obtain the desired effective low energy Lagrangian above. In this 
case there is only one isosinglet (5 = 2/3 leptoquark present, which turns out to be quite 
advantageous. 

In this paper we will consider and classify all models wherein heavy fermions are used 
to generate the effective interactions Cwanted or Chanted lo"^ energies. As we will see, the 
emphasis of our approach is somewhat different than that of Altarelli et al, in that we will 
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keep the new heavy fermions as active ingredients in our models and not treat them as an 

auxihary device to produce the desired couphng structure. In particular, we will assume that 

exotic, vector-like fermions exist and that the desired interactions are induced through their 

couplings to the leptoquark and their mixing with the SM fermions. The mixing between the 

new fermions and those of the SM will be generated by conventional spontaneous symmetry 

breaking (SSB) via the usual Higgs doublet mechanism. It is only through SSB that the 

above effective Lagrangian can be obtained in the fermion mass eigenstate basis from an 

originally gauge invariant interaction. The small size of the effective Yukawa couplings in the 

above Lagrangians, Chanted oi CJ^antedi then directly explained by the same mechanism 

that produces the ordinary-exotic fermion mixing and automatically sets the scale of the 

vector-like fermion masses in the TeV region. We note that the use of vector-like fermions 

in this role is particularly suitable since in their unmixed state they make essentially no 

contribution to the oblique parameters |jT8[|, they are automatically anomaly free, and they 

can have bare mass terms which are SM gauge invariant. (Alternatively, their masses can 

be generated by the vacuum expectation value of a SM singlet Higgs field.) Mixing with the 

SM fermions does not significantly detract from these advantages as we will see below. As is 

by now well-known the leptoquark itself does not significantly contribute to the oblique 

parameters provided it is either an isosinglet, which will be the case realised in all of the 

models below, or in a degenerate multiplet. 

Before discussing the construction of new leptoquark models with vector-like fermions, 

it is interesting to note that HERA will not be able to distinguish between the two scenarios 
described above, even if the relative ej and vj branching fractions are precisely measured. 
The only means of differentiating the models is to either find the other new particles antic- 
ipated in each scheme, or to directly produce the ~ 200 — 220 GeV leptoquarks at a high 
energy e+e~ collider such as the NLC^]. As we will see below, the charge and weak isospin 
of the leptoquark is fixed in the models with vector-like fermions and is independent of the 
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Figure 1: Cross section(left) and associated polarization asymmetry(right) for the produc- 
tion of a pair of 200 GeV leptoquarks at a 500 GeV NLC. The dashed curve is the model 
of Babu, Kolda and March-Russell while the solid line is the prediction of the model with 
vector-like fermions. 



value of Bi. However, in the Babu et al. approach the leptoquark's effective weak isospin is 
highly correlated with the value of B^. Fig.|l] displays a comparison of the leptoquark pair 
production cross section and polarization asymmetry for these two models at a 500 GeV 
NLC. It is clear that unless Bi is very close to 50% the two scenarios will be easily separated 
at the NLC. These results also show that a leptoquark with the quantum numbers antici- 
pated in vector-like fermion models is trivially distinguishable from the more conventional 
BRW leptoquarks by the same analysislH]. 



1.3 Constraints on Leptoquark Coupling Parameters 

As we will find below, in models with vector-like fermions, the only new physics at low 
energies introduced by the leptoquark itself can be parameterized in terms of the interactions 
in Chanted- It is then straightforward to use existing data to constrain the effective Yukawa 
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couplings \u,d', here, we can express A„ in terms of Bi = A^/(A^ + A^), since we assume that 
the leptoquark has no other decay modes. As discussed above, the Tevatron searches place 
a Xd independent constraint on Bi for any fixed value of the leptoquark mass. Similarly, the 
recent measurements of Atomic Parity Violation (APV) in Cesium [0] place Bi independent 



bounds on AdpOf for fixed values of Mlq. In addition, universality in vr decay constrains 
the product AuA^I^, while the observed rate of NC events at HERA constrains instead the 
product X'^Bf, in the later case QCD corrections are quite important [^]. The latest available 
results presented by both the ZEUS and HI Collaborations in the neutral current as well 
as the charged current are included in our estimate of the cross section constraints for both 
channels. (We note that due to the relatively low statistics and other uncertainties the errors 
in this case are probably significantly underestimated so that this band is actually somewhat 
wider than what is shown below.) Combining these constraints defines an approximate 

allowed region in the Bi — A^ plane which is presented in Fig.^ for Mlq = 200, 210, 220 GeV. 
Here, A = A/e with e being the conventional proton charge (this scaling of the coupling to 



e follows earlier tradition p^). We note that the size of the (apart from the HERA data) 
95% CL allowed region is sensitive to the two possible choices of the sign of the product of 

XuXd- As we will see below, the region corresponding to XuXd > is preferred so that the 
vr decay data has little impact in restricting the parameter space. From Fig.|^ we see that 
the position of the allowed region moves up and to the right as the mass of the leptoquark 

increases from 200 to 220 GeV. For the case AuA^ > 0, the size of the allowed region is not 
greatly affected as the leptoquark mass increases whereas, for XuXd < the region grows 
significantly in area with increasing mass. The size of the allowed region in each case would 
be substantially smaller if CDF and DO could combine their results and further constrain 
the value of Bi. 
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Figure 2: Allowed parameter space region in the Bi — plane for a leptoquark with mass 
200 GeV(top left), 210 GeV(top right) or 220 GeV(bottom). The region allowed by Tevatron 
searches is below the horizontal dotted line while that allowed by APV data is to the left of 
the vertical dotted line. The region inside the solid band is required to explain the HERA 
excess in the NC channel at la. The region between the dashed curves corresponds to the 
la range required to explain the apparent excess at HERA in the charge current channel. 
The region above the dash-dotted curve is allowed by tt decay universality: the lower(upper) 
curve corresponds to the case where A^Ad > (<)0. 
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In addition to the constraints shown in FigM further leptoquark couphng information 



can potentially be obtained [p!^ from examining the sum of the squares of the first row of the 
Cabbibo-Kobayashi-Maskawa (CKM) weak mixing matrix, X)j iK^iP- In the SM this sum is, 
of course, unity, but leptoquark exchange in /? decay can yield either an upward or downward 
shift in the extracted value of \Vud\ of 

\Vud\eff - \Vud\true 1-52 X lU |^ j j^g l^o.isj ' 

SO that it would appear experimentally as if a unitarity violation were occurring. Interest- 



ingly, the value of the above sum has recently been discussed by Buras|25], who reports 
J2i iKiiP = 0.9972 ± 0.0013, which is more than 2a below the SM expectation. Clearly, 
if XuXd > 0, leptoquark exchange provides one possible additional contribution which, for 
Am = \d = 0.15 (implying Be = 1/2) and Mlq = 200 GeV, would increase the sum to 
the value 0.9987. If we take this new determination of Vud seriously, then the constraint on 
leptoquark parameters from CKM unitarity can be written in terms of A^ and at the la 
level as (for the case of same sign leptoquark couplings) 



f200 GeV\ f \d\ I -Be 
2.8 ±1.3 = 1.52 \^] \ 4 

which is easily satisfied over most of the allowed parameter space in Fig.^. As we will see 
below, the mixing between the SM and vector-like fermions can also yield an additional small 
positive or negative contribution to {VudUff which can have an effect on the CKM unitarity 
condition in some models. 
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2 Analysis and Construction of New Leptoquark Mod- 
els 



Employing the BRW assumptions (a)-(e) listed above we can construct our new extended 
set of leptoquark models using the following prescription: 

(i) The leptoquark couples a SM fermion multiplet, one of (L, Q, vP, d!^, e"^, {v"^)), where 
Q is the usual left-handed quark doublet, to an exotic vector-like fermion Xj (or X^) in a 
gauge invariant manner. For simplicity, the vector-like fermion is assumed to be an isosinglet 
or isodoublet under SU{2)l and either a singlet, triplet, or anti-triplet with respect to 
SU{2i)c- will denote the new fields with fermion number F > (<)0. 

{a) If Xj(Xf) couples a SM fermion with a given hehcity to the leptoquark, then 
Xf{Xi) couples via H or to the SM fermion of the opposite helicity, where H/H^ are 
conventional doublet Higgs fields. We introduce both H/ H'^ fields as independent degrees of 
freedom to allow for supersymmetrization of the models we construct. 

{in) To obtain the effective Lagrangian in Eqn. (1) we require that terms of the form 
MU'^ + h.c. and £1)^ + h.c. must both appear in the original Lagrangian before spontaneous 
symmetry breaking by the H/H'^ vacuum expectation values (vevs), where one of HilC^) 
and £{V^) must be a SM fermion field. This insures that a Q = — 1(0) lepton will couple to 
Sk Q = —1/3 (-1-2/3) quark to produce an F = leptoquark and that the type of structure in 
(^wanted Can bc obtained after mixing. 

{iv) Bare mass terms for the fields Xi of the form MiXiX^ must be added to the 
original Lagrangian. 

{y) We follow the BRW assumptions (a)-(e) catalogued in the introduction. 
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We note that in the supersymmetric version of these models, the conjugate leptoquark 
field LQ^ must also be present and that it cannot couple directly to any of the SM fermion 
fields, due to gauge invariance, unless it mixes with the leptoquark. This implies, in the 
zero LQ — LQ'^ mixing limit, that the conjugate leptoquark field cannot be produced at 
HERA, and that its production signature at the Tevatron will necessarily be quite different 
than that of the leptoquark and will have thus escaped detection, even though the LQ^ pair 
production cross section is the same as that for leptoquark pairs of the same mass. We will 
briefly discuss the more complex situation which includes this type of mixing below. 

We now begin to classify all possible models which employ SM and vector-like fermion 
mixing to obtain the desired leptoquark couplings. We will take one SM fermion multiplet 
at a time and pair it with a vector-like fermion and a leptoquark. Since there are six SM 
fermion multiplets (allowing for the possibility of u'^) there are naively at most six possible 
models that can be constructed. (As we will see the actual number is somewhat more than 
this since various combinations of these models are feasible.) To demonstrate how these 
construction rules work in practise, we begin by considering the flrst case in detail. Here, 
we couple an exotic fermion, denoted as Xi, to L plus a leptoquark, i.e., LX^ ■ LQ. In this 
case [iii) above requires that Xi be an isodoublet, with member charges of 2/3, —1/3 since 
the leptoquark charge is fixed, as well as an SU{3)c triplet. The BRW assumption (a) then 
dictates that the leptoquark be an isosinglet. We can thus write Xf = ([/°, D^), where the 
superscript denotes the weak eigenstate fields, (ii) and (iv) above then instruct us to add 
the SM gauge invariant terms Xiu'^H + X^d'^H'^ and M^XiX^ to the Lagrangian. Including 
the Yukawa couplings (which we assume are of order unity) these terms, together with the 
gauge interactions of both the leptoquark and the fermion doublet X^, form our new set of 
interactions that are added to the SM. Denoting this as model A, we thus arrive at 

Ca = ^aLXI ■ LQ + ttuXiu^H + aaXid^H'' - MiX^X^ + gauge + h.c. , (5) 
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where ^ gauge'' represents the new gauge interactions of the leptoquark and Xi. We emphasize 
that all of the above Yukawa couplings are assumed to be of order unity. 

When H and receive vevs {v and v'^), the a^^a terms in the above Lagrangian 
induce off-diagonal couplings in both the Q = —1/3 and Q = 2/3 quark mass matrices. 
Neglecting the u- and d-quark masses, these are given in the 'ip^Mt/j^ weak eigenstate basis 

by 





/ 




\ 


\^ a^v -Ml J 


\ 







R 




(6) 



(7) 



Both My can be diagonalized by a bi-unitary transformation which becomes bi-orthogonal 
under the assumption that the elements of ^ are real, resulting in the diagonal mass 
matrices M*^^ = UL{u,d)Mu,dUR{u,dy . Since UL,R{u,d) are simple 2x2 rotations they 
can each be parameterized by a single angle O^^j^. For the case at hand 9^''^ = 0, whereas 

9^"^ ~ ttu^dviv'^) /Ml. Taking the Yukawa couphngs to be of order unity and v, v'^ ~ 100 — 250 
GeV, the size of the mixing is fixed by the scale of Mi. Writing ~ + 9^u in terms of 
the mass eigenstate fields, and similarly for the interaction involving the SM fermions 
and the leptoquark thus becomes 



light 



' i^u + I — — — I ed 



Ml 



V Ml 



LQ + h.c. , 



which is the exact form we desired in Eqn. (1). This naturally leads to a reasonable relative 
branching fraction for the LQ vj decay mode, and gives acceptable values for Xu^d in Eqn. 
(1) for Ml in the 1-5 TeV range. 
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At this point one may note that we have omitted a term in La of the form —M'QX^, 
with M' being a bare mass parameter. Such a term is, of course, gauge invariant and 
should be present in principle but has little influence on the scenario as far as the leptoquark 
interactions are concerned. Of course one can always invent a symmetry to forbid this term 
if so desired as in practice such a term may produce an uncomfortably large mass for the 
SM fermions, induced by mixing, and so additional care is required. However, to keep the 
following discussion as general as possible, such terms will be included in our discussion. 
With M' being the same order as Mi there is essentially no change in our result for the 
right-handed mixing above; we now obtain 6^"^ ~ au4v{v'^)Mi/ [Ml + M'^) ~ au4v{v'^) / M^. 
However, M' induces a non-zero mixing for the left-handed fields, but this does not infiuence 
either the leptoquark or Z boson couplings to the light fermions. There is a new contribution 
in the case of the light fermions' charged current couplings to the W , but it is quite suppressed 
being proportional to A = 1 — cos(^^ — O'l) , with the difference Ol — O'lc^ (^m^^ — (^d'^"'^) / 
being small. Note that while both 9^''^ are large, neither is directly observable and it is the 
difference between the two, which is very small, that is observable. It is also important to 
remember that this mixing angle difference is also proportional to i?| — (1 — BiY, so that 
it is further suppressed for values of approaching 0.5. Even without considering these 
cancellations, we estimate the effect to be very small since the difference in the left-handed 
mixing angles is roughly given by — ~ ~ (0.05)^, implying A < 10~^. In the next 
section we will return to the general question of whether the effects associated with the finite 
size of these mixing angles can lead to observable shifts in SM expectations. 

To proceed with our systematic analysis, we first list the remaining five skeleton mod- 
els that are obtained by simply combining the other SM representations with an appropriate 
vector-like fermion and leptoquark field (note that both models B and F involve the field 
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— ^bQX2 ■ 


LQ + OeXae^//^ + a^X2v''H - M2X2XI 


Cc 


= XcX^u'' 


■ LQ + aiLX^H - M^X^X^ , 




— \j)X/^d'^ 


• LQ + a2LXlH'' - M^X^Xl , 


jCe 


= XeX^c'- 


LQ + asQX^H'^ - M^X^X^ , 




= XfXqv" 


■ LQ + a^QX^H - M^X^X^ , 



(9) 



where the usual 'gauge + h.c.^ terms have been dropped for simpHcity. Note that model B 
is essentially the leptonic equivalent of model A; here, the vector-like fermion field X2 is a 
color singlet, weak isodoublet, i.e., = {N*^, E^), and the leptoquark remains an isosinglet. 
This model requires the neutrino to be a light Dirac field or, at the very least, to appear 
as missing pt in the leptoquark decay process. 

It is important to notice that some of these individual skeleton models do not satisfy all 
of the model building constraints listed above, in particular {iii). However, this requirement 
can be satisfied by taking combinations of the various skeleton Ci above, taking care not to 
violate the BRW condition (d) that the leptoquark couplings remain chiral. The weaknesses 
in models C and D as well as E and F can be overcome by simply pairing them: 

CcD = [XcNu^ + XDEd""] ■ LQ + a^LN^H + a2LE^H'' - MmNN" - MeEE" , 
Cef = [XeDc" + XfUv"] ■ LQ + a^QD^H" + a^QVH - MdDD" - MuUU" , (10) 

where the superscript '0' denoting the weak eigenstate has been dropped for simplicity. 
Both models CD and EF now satisfy all of our model building requirements, however, it 
is important to realize that these two combinations are not the only set of alternatives. In 
this case, the fields U,D,N and E are identified with X6,5,3,4, respectively, and are all weak 
isosinglets with the field notation designating the color and charge information. Note that 
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the individual bare mass terms are present for these fields, and that in both models the 
leptoquark is again an isosinglet with charge 2/3. As in the case of model B, model EF 
requires v'^ to appear as missing pt in the leptoquark decay, otherwise these models are 
excluded by the apparent HERA CC excess and, possibly, by the Tevatron constraints. We 
note that, as in the case of model A, additional gauge invariant mass/mixing terms can be 
added to the Lagrangians of models B, CD and EF. These take the form of —M'j^LX2 for 
model B, -M'^Nu" - M'^Ee'' for model CD and -M'jjDd'' - M[jUu'' for model EF. They 
produce essentially no additional new physics effects at a visible level in the latter two cases 
since none of the SM fermion couplings to the gauge bosons are further altered. However, 
in model B, as was seen for model A, a modification of the SM leptonic CC couphngs to the 
W boson will occur and is proportional to the square of the difference in the left-handed 
mixing angles needed to diagonalize the neutral and charged lepton mass matrices. For 
completeness, the mass matrices for the vector-like and SM fermion sector for each of these 
models are as follows (using the same weak eigenstate basis as above): for model B, 




(11) 



Me 




(12) 



whereas for model CD we find 



M 




(13) 







a2V 



Me 



(14) 



V 



-M' 



E 
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and for model EF we correspondingly obtain 



Ma 




(15) 



(16) 



where in all cases we have allowed for the additional gauge invariant terms discussed above. 
With the primed and unprimed bare mass terms of roughly the same magnitude, the effective 



are expressible in the same form as in Eq. 8 with the appropriate substitutions of masses and 
Yukawa couplings. It is important to remember that in models CD and EF, which involve 
mixings with isosinglet fermions, the roles of the left- and right-handed mixing angles, 9l 
and 6r, are essentially interchanged with respect to those in models A and B where the 
vector-like fermions are in isodoublets. In all cases the relevant mixing angles are of order 
0.05 as obtained in the case of model A, due to the phenomenological constraints imposed 
on the effective Yukawa couphngs by the HERA data. 

Lastly, we note that models A, B, CD and EF are not the only successful ones that 
can be constructed. We can, e.g., take either model A or B and combine it with one of the 
skeleton models C-F; for example, model B could be coalesced with F. In principle, many 
potential hybrid models of this type can be constructed. This observation will be important 
below when we discuss the unification of these models within a GUT framework, as well as 
the phenomenological implications of the SM and vector-like fermion mixing. We note that 
in these more complex models the fermion mixing (s) that generate the SM fermion couplings 
to the leptoquark can arise from multiple sources. Of course, when we attempt to construct 
further hybrid models, we must take care not to violate the assumption that the leptoquark 




^) couplings can be read off directly from these matrices and the above Lagrangians 
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couplings are chiral. Given this very strong constraint, the entire hst of models that can be 
constructed in this fashion are only ten in number: A, B, CD, EF; AC, AD, ACD, BE, BP 
and BEF. We note that models A, CD, AC, AD and ACD produce the effective interaction 
(^wanted, while models B, EF, BE, BF and BEF produce instead Chanted- The models and 
the exotic fermions associated with each of them are catalogued in Table 3. 

3 Implications and Tests 

Some phenomenological implications of these models are examined in this section. The de- 
tailed phenomenology depends on whether or not supersymmetry (SUSY) is also introduced. 
Clearly, the non-SUSY versions are more easily analyzed but both classes of models share 
many common features which we will discuss here. These include new interactions due to 
the mixing between the vector-like and SM fermions as well as from the existence of the 
leptoquark and the vector-like fermions themselves. 

3.1 Direct Production of Vector-Like Fermions 

The production and decay of vector-like fermions has been extensively discussed in the 
literature p4|, particularly in the context of grand unified theories. The mixing in- 
duced between these new fields and the ordinary SM fermions not only modifies the SM 
fermion couplings to the W and Z but also leads to flavor-changing Z interactions involving 
a single SM fermion and a vector-like fermion. This implies that the vector-like fermions 
can be produced in pairs via the usual mechanisms, or singly via mixing. Once produced, 
they can decay through mixing into a SM fermion and a. Z ot W with comparable rates. 
However, unlike most models containing vector-like fermions, it is more likely here that at 
least some of these states will dominantly decay to leptoquarks instead due to the large 
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Model 



A 



CD 



AC 



AD 



ACD 



B 



EF 



BE 



BP 



BEF 



Vector-like Fermions 



U 
D 



L,R 



/ 

V 

U 
D 



U 
D 

U 
D 



L,R 



L,R 



L,R 



E 



L,R 



L,R 



L,R 



Ul,r', Dlr 



N 
E 

N 
E 



R 



L,R 



L,R 



L,R 



N 
E 



L,R 



Table 1: Listing of models and the vector-like fermions which are contained in them. 
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assumed size of the Yukawa couplings. Given the expected large mass of the new fermions 
in these models, they will only be accessible at the LHC (until ^/s=2-10 TeV e"^e~ or 
colliders are constructed). For masses of order 1 TeV and an integrated luminosity of 100 
fb^^, we estimate the yield of color triplet vector-like fermion pairs at the LHC to be of 
order 10^ events, where they are produced by a combination of gg and qq fusion. If the 
W and Z final states produced in the vector-like fermion decays can be triggered on with 
reasonable efficiency this implies that the production of such heavy states should be rela- 
tively straightforward [^. The production and detection of heavy color-singlet states at a 
reasonable rate seems somewhat more problematic [^] due to background issues. 



3.2 Universality Violations Revisited 

Do the vector-like fermions have visible indirect effects at lower energies? We first examine 
whether the vector-like and SM fermion mixing itself induces a sizeable universality violation. 
Here, the cases where the vector-like fermions are isodoublets or isosinglets induce quite 
different effects; recall that in the isodoublet vector-like fermion scenario, 6ji ~ 0.05 and the 
difference 0^ — Of (0.05)^, whereas the reverse is true in the case of isosinglet fermions. 
We thus find the following shifts in the CKM element in each of the above principle 

models (to leading order in the mixing angles) 

A : -(ei - eif + {eiei)\ 

B : 0, 

CD : 0, (17) 

EF : -{eif-{eif. 

Clearly the effect is very small in the models with isodoublet quarks (model A), but can be 
sizable in the isosinglet quark case (model EF). In fact, for model EF we see that the effect 
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of mixing is to decrease the value of |Kid|ejj relative to the SM expectation by an amount of 

order 10~^; this is at the level of current sensitivity, as discussed in the previous section, and 
is comparable to the size of the present difference between experiment and the expectations 
of unitarity. 

The other consequence to notice above is the null result in the case of models B 
and CD. In these scenarios, the same mixing that affects nuclear decays also appears in 
the calculation of fi decay and is therefore absorbed into the definition of Gp- However, a 
residual effect from the mixing will remain in the ratio of widths for it ^ eu to n ^ fiu. This 
results in another shift in these models, in addition to that from the leptoquark exchange 
discussed above, from the SM expectation for this ratio by an amount 

CD : -{eir-iei)\ (18) 

which is negative and can be sizeable in the case of model CD for mixing angles of order 
0.05. Experimentally , the value of this universality testing ratio to its SM expectation is 
found to be 0.9966 ± 0.0030; we note that the potential deviation from unity is comparable 
to the expectation in model CD. 

Lastly, we note that in model B a right-handed charged current is generated for the 
electron, which could in principle be observed in /i decay if u'^ appeared as missing energy 
or pt- However, the size of the right-handed amplitude generated through this mixing is far 



too small (by several orders of magnitude) to be detected in the Michel spectrum pT 



3.3 g — 2 of the Electron and Electron Neutrino 



One reason for demanding that the leptoquark couplings to fermions be chiral is to avoid 
the enhancement of a number of loop-order processes, e.g., the g — 2 of the electron. Here 
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we have successfully constructed cliirally coupled leptoquark models and hence their contri- 
bution to the electron's g — 2 is very small, however, there remains the possibility that the 
mixings between the SM and vector-like fermions may reinstate significant contributions to 
Og. Model B provides an example of this scenario, since in this case, both left- and right- 
handed leptonic couplings to the ly-boson exist and the heavy can participate in Oe as an 



intermediate state. The contribution in this case can be immediately obtained from Ref. |28 



and directly compared with the prediction of the SM||29] and the experimental value. For 
the difference between the latter two we find (with the total uncertainty in the difference 
given in parenthesis) 

^exp _ ^SM ^ _i3_2(27.2) X 10-12 ^ ^19) 
while the additional contribution in the case of model B can be written as 

Aae = (-34346 x 10-^^)F{x) sin(^^ - 61) sin , (20) 

where F{x) is a kinematical function of the mass ratio x = M'^/M^. The large numerical 

size of the prefactor gives some warning that the effect might be of a reasonable magnitude 

even though it is highly suppressed by several powers of mixing angle factors. Taking 6[i 

~ 0.05 and the difference 01 — 91 ~ (0.05)^ as usual, we obtain the results displayed in Fig.|^; 

note that the absolute value of the shift is presented since the signs of the mixing angles 

are unknown. This analysis demonstrates that for typical ranges of the parameters in this 

model, the size of Aa^ is comparable to or larger than the present uncertainty in the value of 

Qe] hence values of ^ 5 TeV are excluded for these suggestive sizes of the mixing angles. 
In a similar fashion the corresponding contribution to the magnetic moment of the 

electron neutrino, be obtained, provided that Ue is a Dirac fermion. (We recall that 

both the electric and magnetic dipole moments of a Majorana neutrino vanish identically.) 

In this case the amplitude arises from a penguin diagram with the vector-like fermion E in 

the intermediate state. The results thus take similar form to that for Aae above, except that 
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Figure 3: Contribution to the anomalous magnetic moment of the electron in model B in 
units of 10~^^ due to mixing between the SM and the vector-hke fermions as a function of 
the N fermion mass. 
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the kinematic function is different and with the replacement sinOfj^ sin^^. The result of 
this calculation is presented in Fig.^ and should be compared to the present experimental 
limit of < 180 x lO'^'^fi^ at 90% CL from elastic elastic scattering using reactor 
neutrinos. Stronger bounds (by factors of order 10) based on astrophysical constraints remain 



somewhat controversial ||2 1|| . Note that a similar graph without the attached photon is capable 
of generating a mass for z/g in the range 10^'^ — 10~^ eV. 
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Figure 4: Mixing induced contribution to the magnetic moment of the electron neutrino in 
model B in units of 10~^^ Bohr magnetons as a function of the E fermion mass. 
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3.4 Oblique Parameters, Z Pole Observables, and APV 



Vector-like fermions are known to have negligible contributions to the oblique parameters [T 
However, once vector-like fermions mix with their SM counterparts it is possible to induce 
non-zero shifts in the values of these parameters. As a numerical example, we examine the 
size of these contributions in model A. In the case of the shift in the p parameter, Ap = aT, 
there are two sources which contribute here: (i) the modification of both the vector-like and 
SM fermion couplings to the W and Z due to mixing; (ii) the U and D masses, originally 
degenerate, are now split by an amount — M^, = a^v"^ — a^v^"^. Writing in this case 
Cu,d = COS 6^'^ and Su,d = sinO^'^, one obtains 

Note that when c„ = q = 1, Mu = Mo and Ap vanishes as expected. For Mu,d — 5 TeV 
and 6^'^ ~ 0.05, Ap is found to be < 10~^, far too small to be observed. In a similar vein, 
the induced value for the parameter S is found to be less than 5 x 10~^ and is hence also 
vanishingly small. Thus, although the vector-like fermions do not remain purely vector-like 
after mixing, their contribution to the oblique parameters remain negligible. This same 
pattern is repeated in the case of the other models with only minor differences, e.g., color 
factors are present in the case of model B and the gauge invariant mass terms for the two 
isosinglet fields in either models CD or EF can be different. Numerically, however, similarly 
small results are obtained for the oblique parameters in these remaining models. 

Are there observable modifications in the SM fermion couplings to the Z-boson? 
Recall that, for example, the mixing of the u and d quarks with vector-like fermions which 
have weak isospin Tg^g^ produces a shift in the u and d couplings to the Z of Avu{au) = 
{Tl - l/2){sl)' ± r^„(s^)2 and At;,(a,) = (T^, + l/2Mf ± T^,(4)', respectively, using 
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the notation above. (The corresponding shifts in the case of leptonic mixing can be obtained 
from these expressions by trivial notational changes.) Two places where these coupling shifts 
may show up most clearly are in the partial widths of the Z-boson and in APV. In both 
these observables there is an additional shift in the case of the leptonic couplings due to the 
overall change in the coupling normalization from the redefinition oi Gp from muon decay, 
as discussed above. However, the Z leptonic asymmetries, which are particularly important 
observables, are insensitive to these overall changes in the coupling normalization. For this 
case, taking the relevant mixing angle to be 0.05 as usual, we find that the Z partial width 
to the e"^e~ final state is decreased (increased) by an amount of order ~ 0.2 MeV for the 
isodoublet(isosinglet) model. Correspondingly, the apparent shift in the value of sin^^^-'^-^ 
from the asymmetries increases(decreases) by an amount of order ^ 0.0006 for these same two 
cases. Both of these shifts are essentially at the boundary of the current level of sensitivity 
for LEP/SLD measurements pO| . Similarly there is a corresponding shift in the number of 



neutrinos extracted from the measurement of the Z invisible width by ^ 0.005. 

These shifts in the SM fermion couplings can modify the expectations for APV as 
well since the effective weak charge, Qw, directly probes the two products agVu and aeVd in 
addition to the shift in the overall normalization that occurs with leptonic mixing. In the 
case where the SM fermions mix with their leptonic vector-like counterparts, the shift in 
is directly given by 

AQ./Q. = 5p- 2(T^, + l/2)(sl)2 + 2X^,(4)^ , (22) 

where 6p represents the change in the overall coupling normalization and is given to leading 
order in the mixing angles by 

B : {01 - 9lY - {O^nOnf , 
CD : {9iy + {9lf. (23) 
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We find that this fractional shift in is at the 10~^ level for either isosinglet or isodoublet 
leptonic vector-like fermions and is hence clearly too small to be observed. The modification 
could be potentially larger when mixing occurs in the u and d couplings and where there is 
no overall change in the normalization. However we find that the individual contributions 
of the u and d quarks tend to cancel each other instead of adding coherently, leaving, again, 
a relative shift in Qyj at the 10"'^ level. 



3.5 Drell-Yan Production in the e^Ue Channel 

What future constraints can be placed on the leptoquark couplings? We know from earlier 
workf^ that the coupling can be probed in high precision measurements at LEP II in 
e^e~ qq and also at the Tevatron via NC Drell-Yan production. Can future colliders 
also probe the coupling? One possibility is to examine the corresponding CC Drell-Yan 
process at hadron colliders, p^p^ e^z/. In addition to the usual SM VT-boson exchange, 
leptoquarks can also contribute to this process via t-channel exchange involving both the A^ 
and Xu couplings. The subprocess cross section for this reaction is found to be 



da{ud — ^ e Ue) GpM^ 



w 



dz 12tts 



(24) 



where x = GpM^ /2^/2 Tca and z = cos 6*, the parton center of mass scattering angle between 
the incoming quark and the outgoing negatively charged electron; as usual i = — s(l — z)/2 
and u = —s{l + z)/2. Note that there is no interference between the 1^-boson and leptoquark 
exchanges which will make the leptoquark contribution somewhat more difficult to observe 
although the two distributions peak in opposite angular regions. 

There are two useful observables in this case. First, one can examine the transverse 
mass (Mt) distribution beyond the Jacobian peak associated with VT-boson production. For 

26 




200 400 600 800 



Mt (GeV) 




Figure 5: (a) The electron plus neutrino transverse mass distribution assuming jr^^j < 2.5 
and (b) the folded lepton charge asymmetry in the charged current Drell-Yan production 
channel at the 2 TeV Tevatron for the SM (solid curves) and with 200 GeV scalar leptoquark 

exchange assuming A^A^ = 1 (dashed curves). In (b). from top to bottom in the center of the 
figure, the SM curves correspond to Mt bins of 50-100, 100-200, 200-400 and > 400 GeV, 
respectively. Note that for Mt in the 50-100 GeV range there is no distinction between the 
SM result and that with a leptoquark. 
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large values of Mt one would expect an increase in da/dMr due to the leptoquark exchange. 
A second possibility is to examine the leptonic charge asymmetry, A{rii), for the case of 
electrons in the final state as a function of their rapidity. Here A{ri£) is defined as 

^ _ dN+/d7]i - dN_/dr]e 
~ dN^/dvi + dN./dve ' ^ ' 

where N± are the number of positively/negatively charged electrons of a given rapidity. In 
the SM, the charge asymmetry is sensitive to the ratio of u-quark to d-quark parton densities 
and the v — a structure of the W decay Since the decay structure of the W has been well- 



measured elsewhere 1 32], any observed deviations from SM expectations in this asymmetry 



have been attributed to modifications in the parton density functions ||33||. The possibility of 
new physics contributing to this channel has been overlooked. In calculating the asymmetry 
it is essential to split the integration over the parton densities into 2 regions, corresponding 
to positive and negative lepton rapidities in the W center of mass frame, according to the 
prescription in Ref. [Q. Fig.^ shows how both the binned transverse mass distribution 
and the lepton charge asymmetry, for four bins corresponding to 50 < Mj- < 100 GeV, 
100 < Mt < 200 GeV, 200 < Mt < 400 GeV, and 400 < Mt GeV, are modified by the 
presence of a 200 GeV leptoquark with, for purposes of demonstration, A„Ad = 1. We see 
that the transverse mass distribution does rise above the SM expectations for large values of 
Mt as expected, and that the lepton charge asymmetry can also be significantly modified for 
larger values of Mt- Note that there is little deviation in the asymmetry in the transverse 
mass bin associated with the W peak, 50 < Mt < 100 GeV, so that this Mt region can still 
be used for determination of the quark densities. 

We now perform a analysis to determine the potential sensitivity to leptoquark ex- 
change at the main injector. As shown in Fig.^(a), we divide the transverse mass distribution 
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into several bins corresponding to 



7 bins in steps of 5 GeV in the ranj 


^e 


50 < Mr < 85 


GeV 


10 bins in steps of 20 GeV in the ranj 




85 < Mt < 285 


GeV 


6 bins in steps of 40 GeV in the ranj 




285 <Mt< 525 


GeV 


2 bins in steps of 100 GeV in the ranj 


le 


525 <Mt< 725 


GeV 


1 bin for the ranj 


le 


725 < Mt 


GeV 



This ensures that adequate statistics are maintained in each bin. The apparent rise in the 
cross section in Fig. ^(a) at Mt = 285 GeV is an artifact of the increased bin width at that 
point. For the lepton charge asymmetry the lepton's rapidity is binned as 



12 bins in steps of Arji 


= 0.2 
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range 


50 < Mt < 100 


GeV, 


12 bins in steps of Arji 


= 0.2 


in 
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GeV, 


6 bins in steps of Arji 


= 0.3 


in 
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200 < Mt < 400 


GeV, 


2 bins of \Ar]i\ < 0.4 and \Ar]e\ 


> 0.4 


in 


the 


range 


400 < Mt 


GeV, 



subject, of course, to the constraint Mt < v^e~'^^'. We note again, that there is no sensitivity 
to the leptoquark exchange on the W transverse mass peak (50 < Mt < 100 GeV), and hence 
this region does not contribute to the distribution. The bin integrated cross section and 
asymmetry are then obtained for the SM and for the case of 200 GeV scalar leptoquark 
exchange. We sum over both e+i^e and e~Ue production for the cross section and employ an 
electron identification efficiency of 0.75. The statistical errors are evaluated as 6N = \/N 
and 5A = — A'^)/N, as usual. The resulting 95% C.L. bound in the Bi — plane is 

presented in Fig. ^ for 2 and 30 fb~^ of integrated luminosity with ^/s = 2 TeV. We see that 
even for 30 fb~^, the constraints are inferior to those obtained from present data on tt decay 
universality as shown in Fig. 2. 
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Figure 6: 95% CL bound on as a function of from a fit to both the Mr distribution and 
A{r]e) at the 2 TeV Tevatron for two integrated luminosities as indicated. The area below 
and to the right of the curves are excluded. 
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3.6 Like-Sign Leptoquark Production at the Tevatron 



In models B and C where the u quark couples to a heavy neutral vector-like fermion, A^, 
new processes may arise if is a Majorana field. (Note that for simplicity we have only 
considered the Dirac case in the above discussions.) One such unusual possibility is the 
production of pairs of identical leptoquarks in hadronic collisions via u— or t— channel 
exchange which generates the process uu 2LQ. The leptoquarks then decay to like-sign 
charged leptons plus jets, a relatively clean signature at a hadron collider. Recall that the 
relevant Yukawa coupling involved in this AL = 2 reaction is of order unity so that this cross 
section may be significant even though it is a valence times sea-quark density process at the 
Tevatron. We find the subprocess cross section to be 



where A ~ 1, 2; is defined in the previous section, and Mat is the mass of the neutral vector- 
like fermion. Here, i = — s(l — Pz)/2 and u = — s(l + Pz)/2, where /3 = (1 — Am^Q/sY^'^. 
Note that as Mn —>■ the cross section vanishes as expected for a Majorana fermion induced 
process. The rate for this reaction at the Tevatron with ^/s=2 TeV is shown in Fig.^ for 
A = 1. Here we see that the cross section initially rises with increasing M^v but then begins 
to fall, scaling like as Mjy 00 ■ For Mn=1 TeV, this cross section corresponds to ^ 
100 events at the Main Injector before leptoquark branching fractions are taken into account. 
At the y^=1.8 TeV collider the cross section is smaller by a factor of ~ 0.56. 

3.7 Speculations on a Realistic Flavor Coupling Structure 

Although one can impose discrete or other symmetries so that leptoquarks only couple to 
a single generation in the weak eigenstate basis it is difficult to understand how this might 



da 

dz 



1287r 



Mjv(t + M-2M^) 
{i-M%){u-Ml) 



(28) 
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Figure 7: Cross section for like-sign pairs of 200 GeV leptoquarks at the 2 TeV Tevatron 
as a function of the neutral vector- like fermion mass, Mjv- The Yukawa coupling is assumed 
to be unity. 
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hold in the physical basis. This issue is a major stumbling block for the construction of 
realistic leptoquark models and is one that we have carefully avoided until now. Of course 
the detailed exploration of possible solutions to this problem lies outside the scope of this 
paper[35], however, there are directions that do show some promise |36[ . 



To be more specific, let us concentrate on models which yield the interaction Cwanted 
in Eq. (1) where the u'^ field is absent, and also neglect the possibility of any large leptonic 
mixing. This implies that all of the traditional flavor changing neutral current (FCNC) 
constraints are only to be applied to the quark sector of the model as lepton number is 
conserved. Interestingly, in this case the relevant flavor changing terms are induced by the 
right-handed unitary matrices U{u,d)R of which there is no information since they play 
no role in SM interactions. For purposes of demonstration we will this simply assume that 
element for element they are numerically similar to the corresponding CKM matrix elements. 
Thus flavor mixing now leads us to make the substitutions ur J2i[U{u)R]ii{uR)i and 
similarly for dn in Cwanted- This particular form guarantees that tree level s ^ d and b ^ d,s 
transitions will be accompanied by e~^e~, while c ^ u processes are accompanied only by 
Uei^e- Thus leptoquarks will not mediate the potentially dangerous processes K -kvv or 
D 7re+e~ at tree level. 

Are the tree-level rates induced by these leptoquarks dangerously large? Fortunately, 
the chirality of the leptoquark couplings automatically reduces the size of their potential 
contributions to rare processes and the fundamental couplings present in the Lagrangian are 
already quite small. As we saw from our discussion of \Vud\i for example, the effective Fermi 
coupling for leptoquark exchange was below the level of 1Q~^Gf for typical values of the 
Yukawa couplings. In fact, using A^^^ ~ 0.15 and M^q = 200 GeV it is easy to show that 
this class of models indeed satisfies all of the FCNC experimental constraints in Ref. 



for values of [U{d)R]sd, [U{u)ji]uc of order 0.1-0.2. (This result remains true even when these 
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constraints are strengthened by more recent experimental results |]2T[]). This observation leads 
us to believe that leptoquarks of the type under discussion here are not only compatible with 
present bounds from flavor changing data but may lead to new effects in flavor physics that 
are comparable in magnitude to SM contributions and can thus be searched for in charm or 
B factories Is! 



4 Unification? Never Break the Chain 

4.1 Non-SUSY Case 

If leptoquarks exist and we also believe that there is experimental evidence for coupling 
constant unification then we must begin to examine schemes which contain both ingredients 
as pointed out in our earlier work|^. In the scenarios at hand, the SM quantum numbers of 
the leptoquark are fixed but new vector-like fermions have now been introduced as well, all 
of which will alter the usual RGE analysis of the running couplings. 

Before discussing supersymmetric models we note that coupling constant unification 
can occur in leptoquark models containing exotic fermions even if SUSY is not introduced 



as was shown many years ago in Refs. [37, 38|. Of course in the work of Murayama and 



YanagidaP?!, the leptoquark was an isodoublet which was one of the BRW models, and is 
now excluded by the combined HERA and Tevatron data. In the scenarios presented here 
the leptoquark is now a Q = 2/3 isosinglet so that the Murayama and Yanagida analysis 



does not apply. Fortunately, we see from the results of Ref. that a second unification 
possibility does exist for just these types of models: in addition to the SM spectrum, one 
adds a leptoquark and its conjugate as well as a vector-like pair of color-triplet, isodoublets 
together with the field H'^. This is just the particle content of model A. To verify and update 
this earlier analysis, we assume for simplicity that all the new matter fields are introduced at 
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the weak scale and take sin^ 0^ = 0.2315 as input to a two-loop RGE study. The results are 
shown in Fig.^ where we obtain the predictions that coupling unification occurs at 3.5 x 10^^ 
GeV and as{Mz) is predicted to be 0.118. If unification does indeed occur we can estimate 
the proton lifetime |^9| to be Tp = 1.6 x 10^^^^ years, safely above current constraints pi 
We find this situation to be intriguing and we leave it to the reader to ponder further. 
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Figure 8: Two- loop RGE evolution of the model with the SM particle content together with 
a leptoquark and its conjugate as well as with the vector-like fermions and Higgs content of 
model A. The SU{3)c{SU{2)l, f/(l)y) coupling corresponds to the dotted(ciashed, solid) 



curve. 
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4.2 SUSY Models 



Of course there are other reasons to introduce SUSY beyond that of couphng constant 
unification, so we now turn to the SUSY versions of the above leptoquark models with the 
assumption that i?-parity is preserved, i.e., the HERA excess is due to a leptoquark and not 
a squark produced through R-parity violating interactions. This subject was discussed at 
some length in our earlier work|p from a somewhat different viewpoint but from which we 
are reminded of several important points: 

(i) To trivially preserve the successful unification of the SUSY-SM, only complete 
SU{5) representations can be added to the conventional MSSM spectrum. As is well-known, 
the addition of extra matter superfields in complete SU{5) representations delays unification 
and brings the GUT scale much closer to the string scale. Of course, there still remains the 
rather unnatural possibility of adding incomplete, but 'wisely chosen', split representations. 
This is what happens, of course, in the case of the usual Higgs doublets and is the basis 
for the famous doublet-triplet splitting problem. Employing split representations certainly 
allows more flexibility at the price of naturalness but still requires one to choose sets of 
SU{3)c X SU{2)l X U{1)y representations which will maintain asymptotic freedom and 
perturbative unification. Of course one would still need to eventually explain why these 
multiplets were split. An example of this rather bizarre scenario is the possibility of adding 
a (2,3)(l/6) from a 15 and a (1,1)(1) © (1,3)(— 2/3) from a 10 to the spectrum at low 
energypj. Here the notation refers to the {SU{?>)c, SU{2)i){Y /2) quantum numbers of the 
representation. We remind the reader that in this notation the leptoquark itself transforms 
as (l,3)(2/3); the smallest standard SU{J)) representation into which the LQ + LQ^ can be 
embedded is a 10©10, while in flipped-S'f/(5) x f/(l)01, it can be placed in a 5©5. 

{ii) Since we only have vector-like fermions in our models, it is clear that only pairs of 
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representations, R+R, can be added to the MSSM spectrum in order to maintain anomaly 
cancelation. Of course this is also true for the leptoquark superfield in that both the LQ 
and LQ*^ fields must now be present as discussed above. 

[iii) To preserve perturbation theory and asymptotic freedom up to the GUT scale 
when adding complete representations, at most one 10+10 pair or three 5+5 pairs can be 
appended to the low energy spectrum of the MSSM apart from SM singlets. The reason 
for this is the general observation that if one adds more than three, vector-like, color triplet 
superfields to the MSSM particle content then the one-loop QCD beta function changes 
sign. Recall that the leptoquark itself already accounts for one of these color triplets. This 
same consideration also excludes the introduction of light exotic fields in higher dimensional 
SU{3)c representations. Complete SU{5) representations larger than 10+10 are found to 
contribute more than this critical amount to the running of the QCD coupling which would 
then blow up long before the GUT scale is reached. Whether unification with strong coupling 



is possible has been considered elsewhere]^, but we disregard this possibihty here 



These are highly restrictive constraints on the construction of a successful GUT sce- 
nario containing both vector-like fermions and leptoquarks and we see that none of the 
models discussed above can immediately satisfy them unless the leptoquark and vector-like 
fermion superfields can be placed into a single SU (5) representation. In the standard SU{5) 
picture, we can then place {U, D)^, an isosinglet E'^ and LQ'^ into a single 10 with the 
corresponding conjugate fields in the 10. This would form a hybrid of model A with the 
'skeleton' model D, which we have denoted by AD in Table 1. Of course we pay no penalty 
for also including the 'skeleton' model C here as well, which then yields model ADC. Instead, 
when we consider the flipped-S't/(5) x f/(l) case, it would appear that we can place (A^, E)'^ 
and LQ'^ into a 5 with the conjugate fields in the 5; this is exactly model B. It would also 
seem that no penalty is paid as far as unification is concerned for including the 'skeleton' 
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model C here as well except that this would violate our assumption about the chirality of 
leptoquark couplings to fermions. However, this model is no longer truly unified since the 
hypercharge generator is not fully contained within the SU (5) group itself and lies partly in 
the additional U{1). While the SU{3)c and SU{2)l couplings will unify, f/(l)y will not join 
them even when arbitrary additional vector-like singlet fields are added. Thus unification 
no longer occurs in this scenario so that this possibility is now excluded. 

The leptoquark embedding situation becomes more perplexing if the leptoquark and 
vector-like fermions cannot occupy the same GUT multiplet. In this case unification and 
asymptotic freedom constraints become particularly tight and we are forced to consider the 
split multiplet approach mentioned above. This means that we add the fields (2, 3) (1/6) © 
(1, 1)(1) © (1,3)(— 2/3) and their conjugates at low energies but constrain them to be from 
different SU{^) representations. In this case the combination (l,3)(2/3) © (1,3)(— 2/3) 
corresponds to the isosinglet leptoquark and its conjugate so what remains can only be the 
vector-like fermion fields. Note that we have again arrived back at models AD and ADC. Are 
these the only solutions? We have performed a systematic scan over a very large set of vector- 
like fermions with various electroweak quantum numbers under the assumption that they are 
either color singlets or triplets, demanding only that {€) QCD remains asymptotically free 



and {a) the model passes the so-called "B-test" which is highly non-trivial to arrange. 
Essentially the B-test takes advantage of the observation that if we know the couplings at 
the weak scale and we demand that unification takes place somewhere then the values of 
the one-loop beta functions must be related. Note that it is a necessary but not sufficient 
test on our choice of models but is very useful at chopping away a large region of parameter 



space. Using the latest experimental data[30|, we find that 



B=^-^ — ^ = 0.720 ±0.030, (29) 
62 - hi 
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where the ±0.030 is an estimate of the corrections due to higher order as well as threshold 
effects and the 6, are the one-loop beta functions of the three SM gauge groups. Note that 
Bmssm = 5/7 ~ 0.714 clearly satisfies the test. If we require that (i) and (u) be satisfied 
and also require that the unification scale not be too low then only the solutions described 
above survive after examining > 7 x 10^ combinations of matter representations. While 
not completely exhaustive, this search indicates the solutions above are fairly unique. It is 
interesting to observe that models constructed around model A produce successful grand 
unification both with and without SUSY. 

Finally we need to briefly comment on the possible relationship between the LQ and 
LQ'^ masses and their SUSY partners. In these SUSY models one might imagine that the 
fermionic partner of the leptoquark, the leptoquarkino, may have a mass comparable to the 
vector-like fermions, i.e., of order 1-5 TeV or so. Why then is the leptoquark itself so light? 
One possible mechanism, discussed in another context by Deshpande and Dutta[0|, is to 
envision a large mixing between the leptoquark and its conjugate that produces a see-saw 



effect analogous to what happens in light stop quark scenarios ||43|| . This possibility will not 
be pursued further here. 



5 Summary and Conclusions 

In this paper, we have obtained a general framework for the construction of new F = 
scalar leptoquark models which go beyond the original classification by Buchmiiller, Riickl 
and Wyler. This approach is based on the observation that in any realistic extension of the 
SM containing leptoquarks it is expected that the leptoquarks themselves will not be the 
only new ingredient. This construction technique is, of course, far more general than that 
required to address the specific issue of the HERA excess and, as outlined, can also be used 
to obtain a new class of F = 2 scalar leptoquark models if so desired. 
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To extend leptoquark models into new territories it is necessary to re-examine the 
assumptions tiiat iiave gone into the classic BRW framework. While the assumptions of 
gauge invariance and renormalizability are unquestionable requirements of model building, it 
is possible that the other conditions one usually imposes are much too strong-unless they are 
clearly demanded by data. This observation implies that for leptoquarks to be experimentally 
accessible now, or anytime soon, their couplings to SM fermions must be essentially chiral and 
separately conserve both Baryon and Lepton numbers. The assumption that leptoquarks 
couple to only a single SM generation is surely convenient by way of avoiding numerous 
low energy flavor changing neutral current constraints but is far from natural in the mass 
eigenstate basis. Our analysis indicates that the natural imposition of this condition in 
the original weak basis, and then allowing for CKM-like intergenerational mixing does not 
obviously cause any difficulties with experimental constraints, especially if lepton generation 
number is at least approximately conserved. What is required to obtain a new class of 
leptoquark models is that the leptoquarks themselves must be free to couple to more than 
just the SM fermions and gauge fields. 

Given the fixed gauge structure of the SM the most hkely new interactions that 
leptoquarks may possess are with the Higgs field(s) responsible for spontaneous symmetry 
breaking and with new vector-like fermions that are a common feature in many extensions 
of the SM. Such particles have the advantages that are automatically anomaly free and give 
essentially no significant contributions to the oblique parameters. In the analysis presented 
above we have shown how two particular new forms of the effective interactions of leptoquarks 
with the SM fermions, consistent with Tevatron searches, the HERA excess in both the NC 
and CC channels and low-energy data, can arise through the action of vector-like fermions 
and ordinary symmetry breaking. The typical vector-like fermion mass was found to lie in 
the low TeV region and they could thus be directly produced at future colliders with known 
rates. With our set of assumptions, we obtained ten new models which fell into two broad 
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classes according to the chirality of the resulting leptoquark couplings to the SM fermions. 
The vector-like fermions themselves were shown to lead to a number of model-dependent 
effects which are close to the boundary of present experimental sensitivity including (i) 
violations of quark-lepton universality (for which, as discussed, there is some evidence at the 
2(7 level arising from the CKM matrix), (ii) possible small changes in the Z-pole observables 
for electrons, (in) a small contribution to the shift in the value of the weak charge measured 
by atomic parity violation experiments over and above that induced by the leptoquark itself, 
(iv) a new contribution to the anomalous magnetic moments of the electron and electron 
neutrino, and (v) the possible production of like-sign leptoquarks with a reasonable cross 
section at the Main Injector. We also showed that, as in the case of Drell-Yan in the e+e" 
channel at the Tevatron discussed in our earlier work, there is some potential sensitivity 
to channel leptoquark exchange in the corresponding e^^u channel through the transverse 
mass distribution and the charged lepton asymmetry. 

Leptoquarks within the framework of models containing vector-like fermions were 
shown to be consistent with Grand Unification in both a supersymmetric and non-supersymmetric 
context. The common feature of both schemes is the structure associated with model A, i.e., 
the vector-hke fermions are color triplet, weak isodoublets in a (2, 3) (1/6) representation and 
both H and H'^ Higgs fields are required to be present as is LQ'^ field. In both scenarios 
the GUT scale is raised appreciably from the corresponding model wherein leptoquarks and 
vector-hke fermions are absent. In the SUSY case a (1, 1)(1) field is also required with the 
optional addition of a SM singlet, corresponding to models AD and ACD. In some sense, 
ACD is the "anti-E^g" model in that the color triplet vector-like fermions are in isodoublets 
while the color singlet fields are all isosinglets. Interestingly, in this scenario there is a 
vector-like fermion corresponding to every type of SM fermion. 

Realistic leptoquark models provide a rich source of new physics beyond the Standard 
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Model. 
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